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Abstract 

A detailed derivation of a two dimensional (2D) low energy effective model for 
spinless fermions on a square lattice with local interactions is given. This derivation 
utilizes a particular continuum limit that is justified by physical arguments. It is 
shown that the effective model thus obtained can be treated by exact bosonization 
methods. It is also discussed how this effective model can be used to obtain physical 
information about the corresponding lattice fermion system. 



> ' 1 Introduction 

m 
in 

q ■ One powerful approach to ID lattice fermion systems is to utilize a particular continuum 
limit to derive a low-energy effective model that can be solved analytically. For the simplest 
case of spinless fermions with local interactions away from half-filling, this approach leads 
to the so-called Luttinger model [1] which can be solved exactly using bosonization [2]; 
see [31 IU E] for closely related pioneering papers. Exact solubility means a lot in this 
case: not only the partition function but all Green's functions of the Luttinger model can 
be computed exactly by analytical methods; see e.g. [6] [7] and references therein. This 
method can be generalized to ID Hubbard-like systems and is the basis of a paradigm for 
ID interacting fermion systems [8j (for a textbook presentation see e.g. [9J, Sections 29-32). 
In a previous paper we outlined a similar approach to 2D lattice fermion systems [TU]. The 
aim of this paper is to present details and generalization^ of this proposal. Applications 
and additional results are given elsewhere [TTJ [12] . 

Our starting point is the so-called 2D t-t'-V model describing spinless fermions on a 
square lattice with nearest-neighbor and next-nearest-neighbor hopping constants t and t', 
respectively, and with nearest-neighbor density-density interaction strength V/2; see (1201) - 
(|2T|) for a precise definition. We propose a particular continuum limit similar to one allowing 
to derive the Luttinger model from the ID t-t'-V model. Using that we derive a 2D fermion 
model that is a natural 2D analogue of the Luttinger model not only by its relation to 2D 



1 The discussion in [10] was restricted to the special parameter values t' = 0, K = 1/2 and v a = 1/2 
(these parameters are explained in the next paragraph and after ([5]) and ([9]), respectively). 
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lattice fermions but also because it can be treated by exact bosonization methods. However, 
different from the ID case, we find that only parts of the fermion degrees of freedom of 
this 2D Luttinger mode{§ can be bosonized and thus treated exactly. We suggest to treat 
the other degrees of freedom by mean field theory, and we argue that there is a significant 
doping regime where this treatment is adequate. We refer to the fermion degrees of freedom 
that can be bosonized as nodal and the others as antinodal. 

The 2D t-t'-V model is a spinless variant of the 2D Hubbard model which is often re- 
garded as prototype model for high temperature superconductors (HTSC) [H]. As discussed 
in Section I3TT] our approach and terminology is partly motivated by the so-called pseudogap 
in HTSC [15] . Another motivation are mean field theory results about 2D Hubbard-like 
models; see [16] and references therein. Our approach was also inspired by renormalization 
group studies of the 2D Hubbard model that predict a truncation of the Fermi surface to 
four disconnected arcs [T71 IB] ; see also [IH] • 

Our work owes much to other ideas on 2D interacting fermion systems discussed ex- 
tensively in the literature, mainly motivated by the HTSC problem. The idea that 2D 
interacting fermions can be understood by bosonization in each direction of the Fermi sur- 
face [20] was advocated by Anderson [13] . One realization of this idea close to ours appeared 
already earlier: Mattis [21] pointed out an exactly solvable 2D model of interacting fermions 
similar to our 2D Luttinger model but without the antinodal fermions; see [22] for other 
work on Mattis' model. Various other implementations of this idea have appeared in the 
literature since then (see e.g. [23J and references therein), but, to our knowledge, they all 
differ in detail from ours. One approach similar to ours is by Luther [24J who proposed 
an effective model for 2D lattice fermions taking into account the 2D square Fermi surface 
in all detail. In Luther's treatment an exactly solvable bosonized model is obtained only 
if one ignores interactions coupling different ID "chains" in Fourier space. Since we use 
a simplified description of the Fermi surface we can avoid such a truncation. Still, our 
description of the Fermi surface is more detailed than in Mattis' model [2T] since we also 
take into account the antinodal fermions. The importance of the antinodal fermions in 2D 
Hubbard-like systems is well-known from various different approaches to the HTSC prob- 
lem, including Schulz' work on Fermi surface instabilities [25] and the so-called van Hove 
singularity scenarios [2"6] . 

A few remarks on our level of mathematical rigor is in order. As we show, the 2D 
Luttinger model is a quantum field theory, i.e. a quantum model with infinitely many 
degrees of freedom, that is mathematically well-defined. However, its precise definition 
includes short- and long distance regularizations whose specification requires a considerable 
amount of notation. We therefore first present a formal definition of this model in which 
these regularizations are suppressed^ This formal definition has heuristic value, and it 
allows a concise summary of computations and results. Our derivation of the 2D Luttinger 
model from the 2D t-t'-V model determines the details of all regularizations and the model 
parameters. This derivation is done by by exact computations, except for the continuum 

2 The name "2D Luttinger model" here and in the following is short for "2D analogue of the Luttinger 
model" , and we do not intend to suggest by this name that this model necessarily has "Luttinger-liquid" 
behavior [13] . 

3 A similar formal definition of the Thirring-Luttinger model [H [1] is common in the physics literature. 
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limit that corresponds to a sequence of approximations marked as A1-A5 in Sections 15.21 
and 15.41 Physical motivations for these approximations are discussed in Section [3j Our 
treatment of the 2D Luttinger model by bosonization can be made mathematically precise. 
However, to limit the length of the present paper, we only show in Appendix |B] that the 
same mathematical results that allow to construct and bosonize the ID Luttinger model 
rigorously [27] are sufficient to do so also in the 2D case. Our bosonization results are 
otherwise derived by theoretical physics arguments. Mathematically precise formulations 
and proofs of the latter results are given elsewhere |12j . 

Section [2] contains a summary of our results, including a formal definitions of the 2D 
Luttinger model and a formal description of its bosonized form. Section [3] discusses physical 
motivations for our approach^ Section H] fixes our notation and gives mathematically 
precise definitions of the 2D t-t'-V model and the 2D Luttinger model. Section [5] contains 
a detailed derivation of the 2D Luttinger model. Section E] discusses the bosonization of 
the 2D Luttinger model, the exact solution of the bosonized nodal model, and a derivation 
of the effective antinodal model. We conclude with complementary remarks in Section [7J 
Some technical details are deferred to three appendices. 

2 Summary of results 

We emphasize that the Hamiltonians and (anti) commutator relations of boson and fermion 
operators in this section are only formal, and they become well-defined only with particular 
short- and long-distance regularizations made precise in Sections 14.21 and 15.51 In particle 
physics parlance, the Hamiltonian below defines a quantum field theory model with UV 
divergences, and is is well-defined only with a UV cutoff a > (lattice constant) and 
a IR cutoff L < oo (square root of system area). The importance of the former cutoff is 
highlighted by the non-trivial scaling of the model parameters and the boson operators with 
a; see (jSj) and (fit)]) below. This scaling and the short- and long-distance regularizations 
are such that physical results computed from the 2D Luttinger model have a well-defined 
quantum field theory limit L/a — > oo (this is shown for various special cases in the present 
paper and in [H1H2], an d we conjecture that it is always true). 

The 2D Luttinger model is formally given by the Hamiltonian H = H n + H a + H na with 




(1) 



the part involving only on the nodal fermions, 



H a = d 2 x ^(:^x)[rc F 0+0_ + c' F (d 2 + + d 2 ) - //^(x) : Wr.oW J-r,o(x)) (2) 



4 Mathematically inclined readers may wish to skip Sections [2] and [3] at first reading. 
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the part with only antinodal fermions, and 

H na = Jd 2 xg 4 ^2 ^r,o(x)J r ', s (x) (3) 

r,r',s=± 

the mixed part. The colons indicate normal ordering, d± = d/dx± with x± the compo- 
nents of the 2D position coordinates x, ?/4fj(x) are fermion field operators with the usual 
anticommutator relations {if) rs (x), ijr r , a /(y)} = S r yS SjS '8 2 (pc — y) etc., and 

J r , s (x) =:^ s (x)^ s (x): (4) 

are the corresponding fermion densities. The model parameters cp, c' F , Vp, Qj and /i are 
given below. 

Our main result is a derivation of the 2D Luttinger model from the 2D t-t'-V model in 
Section [5] using certain Approximations A1-A5. This derivation determines the regular- 
izations needed to make precise mathematical sense of this model. It also introduces two 
parameters k and Q in the ranges 

0<k<1, tt(1 - k)/2 < Q < tt(1 + k)/2 (5) 

determining the size and location of the nodal Fermi surface (see Section |3T2|) . This deriva- 
tion also fixes the model parameters 

v F = 2 V2 sin(Q) [t + 2t' cos(Q)]a, c F = 2ta 2 , c' F = 2t'a 2 (6) 

Mo = -4t' - [At + V{1 - nf] cos(Q) - [W + 2V(1 - k)(22 - 1)] cos 2 (Q) (7) 

^ = 2^ 2 = 2Va 2 sin 2 (Q), ^ = 34 = 2Va 2 (8) 

and the filling factor 

i/=I + (l-K)(f -1)+ K 2 K-I). (9) 

The parameter i/ a is the filling factor of the anti-nodal fermions and can be computed 
self-consistently [TT] . 



Our second result is to show that the 2D Luttinger model is equivalent to a model of 
2D bosons coupled to the antinodal fermions and with interactions at most quadratic in 
the bosons field operators (Section 16.11) . A mathematically precise statement of this result 
is given in proposition 16.11 It can be formally summarized as follows: the following linear 
combinations of fermion densities 

9A(x) = Jf[J +, a (x) + J _,,(x)] 

nr (10) 



with 



dot v2 , 
a = a (11; 
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define 2D boson operators obeying the canonical commutator relations [II s (x), $ s /(y)] = 
— i5 SjS /<5 2 (x — y) etc. Moreover, the nodal and mixed parts of the 2D Luttinger Hamiltonian 
are equal to 

H n = V fJ d 2 x :(^[(l-7)n s (x) 2 + (l + 7 )(9 s $ s (x)) 2 ] +2 7 9 + $ + (x)9_$_(x)): (12) 



and 



H na = / d 2 x-^= J rj0 (x)<9 s $ s (x 
J yira £ -~'_, 

v r,s=± 



(13) 



with the dimension-less coupling parameter 



1 27r[t + 2fcos(Q)]' V ; 

As we will see, this implies that the 2D Luttinger model is well-defined provided that 
7 < 1. Fortunately this does not restrict the model to weak coupling. We also outline how 
to exactly solve the bosonized nodal model (Section l6.2p . 

Our third result is an effective model for the antinodal fermions obtained by integrating 
out the bosons exactly (Section I6.3p . This effective antinodal model includes an effective 
interaction that has a non-trivial time dependence. We propose to approximate this inter- 
action by an instantaneous one (as will be discussed, this approximation is convenient but 
not essential). We find that the effective antinodal model can then be represented by an 
Hamiltonian as in ( 1371) but with g 3 changed to #3 — g e g with 

Q = ^(l-Kja 2 

sin(Q)7r[t + 2f cos(Q) + ^(l-fi:)sin(Q)] 1 ; 

the renormalization of the antinodal coupling by the nodal fermions. Note that all pa- 
rameters in this effective antinodal model scale like a 2 . Mean field results for the effective 
antinodal model are given in [HQ. Green's functions of the nodal model defined by the 
Hamiltonian in (I3T)]) can be computed exactly 



3 Physical motivations 
3.1 Experiments 

We recall experimental facts about HTSC [15] that motivate our approach and the termi- 
nology we use. 

An important parameter determining the physical behavior of HTSC materials is the 
so-called filling factor v defined as the total fermion number divided by the largest possible 
number of fermions (note that our definition is such that < v < 1). HTSC are typically 
insulators at half-filling [y = 1/2), and they have a conventional and essentially 2D Fermi 
surface in the so-called overdoped regime where v is significantly different from 1/2. In 
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Figure 1: Fermi surface e(k) = /i for the tight-binding dispersion relations in f JTB"]) for i! jt = 
—0.2 and fi/t = —0.672. The full and dashed circles indicate the antinodal respectively 
nodal regions, as explained in the main text. The dashed lines indicate the Fermi surface 
arcs we postulate to exist in our approach. 

the latter regime HTSC behave like conventional metals. Our interest is in the so-called 
underdoped regime where v is close to but different from 1/2. In this regime HTSC have 
exotic properties not fully understood. In particular, the Fermi surface in this regime seems 
to be partially gapped [28J: there still exists a so-called underlying Fermi surface that can 
be well described by a tight-binding dispersion relation [29] (an example is given in figured]). 
The visible part of this Fermi surface consists typically of four line segments called nodal 
arcs (enclosed by dashed circles in figure [I]), and the remaining so-called antinodal parts of 
the Fermi surface (enclosed by full circles in figure [[]) are gapped. The //-dependence of the 
size and location of the nodal arcs has been measured for various HTSC materials; see e.g. 
1 29 1 and references therein. 
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3.2 Theory 



The continuum limit we use to derive the 2D Luttinger model is based on the following 
hypothesis about interacting fermion systems that can be justified by renormalization group 
arguments; see e.g. [30], Chapter 12 and 14 or [31], Chapter 4. 

HI: There exists some underlying Fermi surface dominating the low energy physics of the 
interacting fermion model, and one can modify, ignore or add degrees of freedom far 
away from this Fermi surface without changing the low energy properties of the model 
much. 

As a motivation we first recall the well-known ID case [8]. The non-interacting Fermi 
surface of the ID t-t'-V model consists of two points, and according to Hypothesis HI one 
can perform a continuum limit that amounts to the following approximations: linearizing 
the dispersion relation in the vicinity of the Fermi surface points, modifying interactions of 
degrees of freedom away from these points, and partly removing the short-distance cutoffEl 
One thus obtains a continuum model with two fermion branches, the left- and the right 
movers, representing the fermion degrees of freedom in vicinity of the two Fermi-surface 
points. The approximations described above provide crucial simplifications since they lead 
to a model of interacting fermions that can be mapped exactly to a non-interacting boson 
model [8]. 

In 2D the situation is more involved. However, it is still useful to consider a typical 
non- interacting Fermi surface of the 2D t-t'-V model. The latter is defined by the following 
tight-binding dispersion relation 

e(k) = — 2i[cos(aA; 1 ) + cos(a£;2)] — 4i'cos(aA;i) cos(a/c 2 ) (16) 

with k = (ki,k2) the usual momenta in the Brillouin zone —it /a < kj < ir/a and the 
hopping parameters t > and t' (a > is the lattice constant). As a representative 
example we plot in figure [1] the Fermi surface e(k) = [i for some small negative values of 
t'/t and the Fermi energy \i. 

The behavior of the dispersion relation e(k) is qualitatively different in different regions 
of the Brillouin zone: since (ir/a, 0) and (0,7r/a) are saddle points of e(k), no linear ap- 
proximation of the dispersion relation exists in the so-called antinodal regions close to these 
saddle points. However, in the so-called nodal regions close to points k = (Q/a, ±Q/a) and 
(—Q/a, ^fQ/a) for some parameter Q rj tt/2, one can approximate by a linear dispersion re- 
lation; see (J47J). (The antinodal and nodal regions are indicated by full respectively dashed 
circles in figure [TJ) This suggests that the fermion degrees of freedom in these different 
regions play different roles in the interacting model. We thus rewrite the model by intro- 
ducing field operators xprj labeled by "flavor indices" r = ± and s — 0, ± that represent 

5 The remaining short-distance regularization corresponds to the non- locality of the two-body interaction 
of the Luttinger model. In the limit where this interaction becomes local, the Luttinger model becomes 
equal to the massless Thirring model 3 which is more complicated from a mathematical point of view; see 
e.g. [32] for a discussion of this limit. 
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the fermion degrees of freedom in these different regions: r = ± and s = ± correspond to 
the regions close to the nodal points (rQ/a, srQ/a), and r = ± and s = correspond to 
the regions close to the antinodal points (n/a,0) and (0, tt/o), respectively. Before we can 
use HI above to modify the model we need another hypothesis that is partly motivated by 
experiments on HTSC mentioned in Section I3.ll 

H2: The underlying Fermi surface contains nodal arcs that can be well approximated by 
straight line segments. 

To be more specific, we assume that the Fermi surface contains four points (Q/a, ±Q/a) 
and (—Q/a, TQ/ a ) for some parameter Q w ir/2, and that close to each of these points we 
can approximate the exact dispersion relation by the linearized one. As will be seen, this 
corresponds to a Fermi surface containing four straight line segments as indicated by the 
dashed lines in figure [TJ We introduce a parameter k determining the size of the regions 
where we use the linearized dispersion relations and such that the length of these line 
segments is (1 — K)y/2n/a. We emphasize that we do not make any assumption about the 
underlying Fermi surface in the antinodal region, and whether the antinodal regions are 
gapped or not depends on parameters and can be determined by computations [TTj . 

Our approximations leading to the 2D Luttinger model are similar to the ones in ID 
mentioned above. They are (essentially) justified by Hypotheses HI and H2 above. Simi- 
larly as in ID, due to their linear dispersion relation, the interacting nodal fermions can be 
mapped exactly to non-interacting bosons. 

4 Definitions and notation 

We write "=" to emphasize that an equality is a definition. We write a vector k G R 2 in 
components either as k = (k±, k 2 ), kj G M., or as k = k + e + + fc_e_ with 

k ± d = f (h ± k 2 )/V2, e ± d = f (1, ±1)/V2, 

and |k| = \/k\ + k 2 _ . We use symbols x, y for fermion positions, k, k' for fermion momenta, 
and p for differences of fermion momenta. 

We write u k± < c" short for u k + < c and < c" etc., and u k± ^ c" is short for u k + -ft c 
or k_ c (or both)" etc. 

4.1 2D t-t'-V model 

We consider a diagonal square lattice A with lattice constant a > and (L/a) 2 ^> 1 
lattice sites: The lattice sites are x = (xi,x 2 ) = (ani,ano} with integers rij such that 
—L/2 < x± < L/2, and we assume that L/(2y/2a) G NU The momenta we use are 

6 We this somewhat unconventional large-distance regularization to simplify some technicalities later on. 
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therefore restricted to be in the following sets, 

A* = {k G M 2 : fc ±G f(Z + |)} 
A*= f '{pGM 2 : P± e 2 fZ} 

(we use antiperiodic boundary conditions for the fermions). The Brillouin zone (lattice 
Fourier space) for the fermions is 



BZ = 



k G A* : ——< kj < —, 7 = 1,2). (It 
a a J 



The 2D t-t'-V model is defined on a Fock space with fermion creation- and annihilation 
operators ip^(pc) labeled by lattice vectors x and satisfying the usual canonical anticom- 
mutator relations. We denote as |0) the normalized vacuum state annihilated by all ^(x), 
and our normalization is such that {^(x), ip^ (y)} = S^ y /a 2 . The Hamiltonian of this model 
is 

H wv = H + H' (19) 

with the free part 

h = A-n* - y) - MWaVWCy) (20) 

x,yGA 

and the hopping matrix T(x — y) equal to t/a 2 > for nearest neighbor sites (i.e. if 
|x — y| = a), t' /a 2 for next-nearest neighbor sites (i.e. if |x — y| = \/2a), and zero otherwise; 
jj, G M. is the chemical potential. The interaction part is 

H'=J2 "Vx - y)^ t (x)^(x)^ t (y)^(y) (21) 

x,yeA 

with w(x — y) = V/4 > for nearest neighbor sites x, y and zero otherwise. Note that our 
scaling with a is such that t, t', V and /i all have the dimensions of an energy. We assume 
that \2t'/t\ < 1 (otherwise the qualitative features of the dispersion relation are different 
from what we assume in Section and our treatment does not apply). 

The average number of fermions per site is the filling factor (or filling) defined as follows, 

z/^ f (f) 2 (iV), iV=X«V(x)^(x) (22) 

xGA 

where (■) is the ground (or thermal) state expectation value and the fermion number 
operator. Thus < v < 1, and half-filling corresponds to v — 1/2. We sometimes refer to 
v — 1/2 as doping. 

We recall that the model is invariant under the particle- hole transformation ^(x) -h- 
(_l)Oi+z2)/a^t( x ) equivalent to the following change of parameters 

(t, t', is, fi, V) -> (t, -t', l-u,2V-fji,V). (23) 

The model is also invariant under the following (discrete) rotation and parity transforma- 
tions 

11 : x= (x u x 2 ) (-x 2 ,zi), P:x^-x. (24) 
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Our conventions for Fourier transformation are as follows, 

xeA 

with k G BZ. This allows us to write 

H o= ^(|) 2 [e(k)-^t(k)^(k) (25) 

keBZ 

with the dispersion relation in (fl6l) . The interaction part in Fourier space is 

H '= E (T) 8s ( k i' k 2 ) k 2 ) k4)^ t (k 1 )^(k 2 )^(k 3 )^(k 4 ) (26) 

k^eBZ 

with 

«(ki, k 2 , k 3 , k 4 ) = u(ki - k 2 ) ^ (i) 5 ki-k2+k3-k4,27rn/a (27) 

nGZ 2 

and 

Hp) = -^2"[ cos ( a Pi) + cos(ap 2 )]. (28) 

The physical interpretation of our interaction is that it contains all possible scattering terms 
k 4 — > k 3 and k 2 — > ki weighted by a factor u(ki — k 2 ) and otherwise restricted only by 
overall momentum conservation up to umklapp processes. 

We also use fermion density operators 

p(x) = ^ f (x)^(x) 
with the following conventions for Fourier transformation 

p(p) = fl V(x) e~ ipx = ^(ki)^(k 2 ) 5 k 1+ p,k 2 + 27r n/a 

xeA ki,k 2 eBZ neZ 2 

for p G A* such that \p±\ < n/a. This allows us to write 

H' = J2 ^(x - y)p(x) P (y) = (i) 2 ^(p)p(-p)p(p). 

x,y P 

Our conventions for density operators of the 2D Luttinger model are similar. 

Remark: Note that our normalizations are such that the naive continuum limit a — > 
makes sense, in particular 5 Xjy /a 2 — > <5 2 (x — y) (Dirac delta) and ^2 x a 2 —> f d 2 x (Riemann 
sum). This limit corresponds to the one discussed in Section [71 remark 4. However, we are 
mainly interested in the continuum limit where v is fixed and close to 1/2, and this limit is 
more delicate. Note also that our formulas have a well-defined limit L — > oo, in particular, 
{^(k), ^t( k ')} = [L/(2it)}Hw 5 2 (k - k') and £ k (27r/L) 2 / d 2 k. 
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4.2 2D Luttinger model 



The formal definition of the 2D Luttinger model in Section [2] can be made mathematically 
precise as follows. 

The 2D Luttinger model is an effective model for the 2D t-t'-V model depending on two 
more parameters k and Q in the ranges given in (jSj) and such that 



kG^(N+|), Q e —j — N (29) 



(it would be easy to drop the restrictions in ( l29l) . but they simplify some formulas below, 
and they become irrelevant for L/a 3> 1). 

We introduce fermion operators ipr}(k) labeled by two indices s = 0, ± and r = ± and 
by momenta k in different Fourier spaces A* as follows, 

A* =(ke A* : \k ± + j-\ < 



y/2, 



a 



A* s=± = kGA*: \k. s + l\< 



1 — Km 



V2a 



Note that there are infinitely many nodal (i.e. s — ±) fermion degrees of freedom (since 
there is no restriction in A* =± on k s ), but only finitely many antinodal (i.e. s = 0) ones. 
These fermion operators obey the usual canonical anticommutator relations normalized 
such that 

{^ )S (k),$, )S ,(k')} = V<V<*k,k'(f ) 2 , (31) 

and they are defined on a Fock space with a non-trivial vacuum Q ( "Dirac sea" ) character- 
ized by the following conditions, 

?/V iS=± (k)fi = $ ii=± (-k)n = Vk such that rk s > (32) 

and (Q,Q) = 1, with (•, •) the Fock space inner product. For the antinodal fermions we 
only assume that antinodal filling is the same for r = + and r = — , i.e. 

v * = 7^)-2 E (T) 2 (^^J,o(k)^o(k)^) (33) 



is independent of r = ±. The latter assumption could be easily dropped, but it simplifies 
some formulas, and it can be checked at the end of our computations. Note that < v a < 1, 
i.e. v a defined in ( )33l) is the filling factor for the antinodal fermions. It is not a free parameter 
but can be computed self-consistently [TT]. We denote by colons normal ordering with 
respect to Q as usual, 

:A:=A-(n,A£l) (34) 

for fermion bilinear operators A = ipl s (\£)ip r s (k') etc. (boson operators are normal ordered 
in the same way in Section [6]). We also introduce the cutoff functions 

Jl if \Vs\<^- a and |p_.|<^g= 

X-=±(P) = \ n „ . V ( 35 ) 

10 otherwise. 
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The 2D Luttinger model is defined by the Hamiltonian H = H n + H a + H na with 
^=EE(l) 2 E r ^ :^(k)^ ri .(k): 

s=± keA* r=± 

+ 2 S (i) 2 (^5Z^(p)^(-P)^-' s (p) + ^ 2 S x+(p)x-(p)4+(-p)^-',-(p 

peA* s=± r,r'=± 

(36) 

the nodal part, 

kA* r=± 

^ • (37) 

+ 2^(i)^3J +> o(-p)J-,o(p) 

peA* 

the antinodal part, and 

H na = (i) 2 ^4^4o(-p)^X S (p)^,s(p) (38) 

P eA* r > r ' s=± 

the mixed part. The interactions involve the following Fourier transformed and normal 
ordered fermion densities, 

4o(p) = E(f) 2 :^,o(k-p)^,o(k): 

^ _ (39) 
J r ,s=±(p) = J2(t) 2 J2 :< s (k-p)^, s (k+v^(l-K;)7me_ s /a): . 

keA* n&L 

The model parameters are given in ([6])-(j9]). As mentioned, the 2D Luttinger model is 
well-defined if 7 in (I14p satisfies the condition 7 < 1, and we always assume this is the case. 

As discussed in more detail below, normal ordering is essential to make the terms in- 
volving nodal fermions well-defined, but for the antinodal fermions it only amounts to 
subtracting finite constants. Note that J r ,o(p) is non-zero only for \p±\ < y/2nK/a, and 
J r s=± (p) is defined for p in 

0^={p 6 A- = W <^}. (40, 

The relation between the precise definition of the 2D Luttinger model here and the 
formal one in Section [2] is spelled out in Section 15.51 
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5 Derivation of the 2D Luttinger model 



5.1 Eight flavor model 

We divide our Brillouin zone in eight non-overlapping regions as indicated in figure |2j For 
that we define eight vectors Q r s with r = ± and s = 0, ±, 2 as follows 

Q + , = (7r/a,0), Q_, = (0, vr/a) 

Qr,s = rsQ/a) for r = ±, s = ± (41) 

Q_, 2 = (0,0), Q +i2 = (7r/a,7r/a) 

for some Q ~ 7r/2 obeying the second condition in §5§ Q We also introduce eight corre- 
sponding rectangular regions A* s so that every vector in BZ can be written uniquely as 
Q r>s +k-|-27rn/a for some r G {+, — }, s G {0, ±, 2}, k G A* s , and n G Z 2 (the umklapp term 
2-KYi/a is sometimes needed to translate to the region — 71/ a < kj < 71 /a). Mathematically 
these regions can be defined as follows 

A; o = |k G BZ 

a;, s=± = {k G BZ 

A; 2 = |k G BZ 

for some parameter n satisfying the first condition in ((5])|3 (see figure [2]). Thus our division 
of the Brillouin zone depends on two parameters Q and k where the former specifies the 
locations and the latter the size of the Fermi surface arcs. It is easy to see that, for geometric 
reasons, k and Q are restricted as in Note that A* is the same for r = + and r = — 
and equal to Ag in f l3"U|) . 

We define 

T/V >s (k) = </>(Q r , s + k + 27rn/a) Vk G A; s 

with n G Z 2 such that Q rs + k + 27rn/a G BZ. This allows us to rewrite the free- and 
interaction parts of our Hamiltonian in ( 125]) and (126]) as follows, 

H ° = E E (t) 2 HQ^ + k) - A*]^,.(k)^(k) 

r,s keA* s 

and 

#' = E E (f) 8 ^(^1^2,^3,^4)^, Sl (k 1 )^ 2iS2 (k 2 )^ 3iS3 (k 3 )^ 4iS4 (k 4 ) 

r^.Sj kjeA*. iS . 

7 This condition ensures that Q r s= ± £ A*. 

8 This condition ensures that mxj{\pla) and (1 — n)?r/(y/2a) are half-integers multiples of 2n/L. 



K7T 



fc. + f + r 



(2Q - tt) 



V2a 



< 



KTC 



V2a 



I < 



(1 - k)tt 
V2a 



(42) 
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Figure 2: Division of the Brillouin zone BZ in eight regions Q rs + A* s marked as (r,s) 
for r = ±, s = 0,±, 2; see (JUJ) and f|42l . The eight dots mark the points Q r ,s, and the 
parameters are Q = 0.457T and k = 0.8. Note that kip can be shifted by integer multiples 
of 27r/a so that the union of these regions cover the BZ exactly once. 

where Kj is short for (kj, 7j, Sj) and 

v(K u K 2 ,K 3 ,Kt) =u(Q n , Sl - Q r2 , S2 +kx - k 2 ) ^ (A) 2 

nez^ (43) 

X ^Qn, si — Qr 2 , 42 +Qr3,a3— Qr 4 ,«4+ki— k 3 +k3— k4,2ir/an- 

The fermion number operator can now be written as 

N = E (x) 2 ^( k )^(k) = £ AT riS) AT r , s = £ ) 2 ^ jS (k)Vi rjS (k) (44) 

keBZ r,s keA*_ s 

where A^. |S is the number operator for the (r, s) fermions. 

Note that, up to now, we only rewrote the lattice Hamiltonian without any approxi- 
mation. However, we now can interpret it as model of eight different flavors of fermions 
distinguished by the labels s and r. The degrees of freedom with the flavor index s = 0, 
r = ± are what we call antinodal fermions, and the ones with s — ±, r — ± are the nodal 
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fermions. We call the degrees of freedom with s = 2 for r = — and r = + the in- and the 
out fermions, respectively. 



We have introduced two parameters Q and k for generality, and the effective model 
we derive below depends on them. In principle these parameters can be determined from 
filling and by minimizing the total energy, but for now we leave them arbitrary. We note 
that the choice k = 1/2 and Q = tt/2 is special since then all the "small" Fourier space 
regions A* s are identical and equal to the Brillouin zone of a "large" lattice with sites 
X = ni(l, l)a + 7i2(— 1, nj integers, which contains eight sites of the original lattice per 
elementary cell. However, it is important to note that our eight flavor model is not local 
on this "large" lattice. 

Our fermion flavors (r, s) can be naturally divided into four different classes with differ- 
ent physical behavior: (i) (+,0) and (—,0) (antinodal fermions), (ii) (+,+), (—,+), (+, — ) 
and (— , — ) (nodal fermions), (iii) (—,2) (in- fermions), and (iv) (+,2) (out-fermions) . It is 
interesting to note that this division of fermions in four classes is also implied by symmetry 
considerations: The (discrete) rotation and parity transformations in (|24p mix our fermion 
flavors (r, s) as follows, 



for r = ±. Thus the four fermion classes above transform under different irreducible 
representations of the group generated by TZ and V . 

5.2 Partial continuum limit I 

We indent to modify the short-distance details of the eight flavor model such that it becomes 
amenable to exact, non-perturbative computations. For that we follow a strategy which 
has been successfully used in ID: We assume that there is some underlying Fermi surface 
dominating the low energy physics, and we can modify, ignore or add degrees of freedom 
far away from this Fermi surface (in the latter two cases we need to correct the definition 
of doping, of course); see Section 1X21 for physical arguments motivating this strategy. 

We expand the dispersion relations of the different fermion flavors (r, s) in Taylor series 
in k±, and we denote the lowest order non-trivial terms as e r}S , 



e(Q r , )=4t', e(Q r ,±) = -4[t cos(Q) + t' cos 2 (Q)], e(Q r , 2 ) = 4(rt - t') (46) 



TZ : ( r> o) -> (-r, 0), (r, ±) -> (±r, =f), (r, 2) -> (r, 2) 
V : (r, 0) -> (r, 0), (r, ±) -> (-r, ±), (r, 2) -> (r, 2) 



(45) 



e(Q r>s + k) = e(Q r , s ) + £ r , a (k) + . . . 



with the dots indicating higher order terms. We find 



and 



£r,o(k) 
e r ,±(k) 
e r , 2 (k) 



—rcpk + k_ — c' F (k^ + k^), 
rvpk±, 



(47) 



(-lrc F + c' F )(kl + kt) 
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with the constants in ([6]). We see that the nodal fermions (s = ±) have dispersion relations 
approximately linear and with a constant Fermi velocity vp, but the dispersion relations of 
the antinodal fermions (s = 0) are quadratic. The in- and out fermions (s = 2 with r = — 
and r = +, respectively) have energies far away from the Fermi energy, and we therefore 
expect that they can be ignored. This will be done below in Approximation A3'. 

Our first approximation is in the free part of the Hamiltonian: 

Al: Replace the exact dispersion relations e(Q rjS + k) in H by the truncated Taylor series 
e(Qr, s ) +e r>s (k). 

This approximation is only essential for the nodal fermions [11], and we use it for the 
other fermion flavors only for aesthetic consistency. 

Our second approximation is for the interaction part of the Hamiltonian (again this 
approximation is only essential for terms involving nodal fermions). 

A2: Replace the interaction vertex fijcfy i n H' by 

v(K l ,K 2 ,K 3 ,K A ) = (^) 2 <5 kl _ k2+k 3_ k4i oM(Qr llSl -Qr 2 , S2 ) 

v X ( 48 ) 

A / j Qr 1 ,3 1 -Qr2,s 2 +Qr 3 ,S3-Qr 4 ,3 4 ,27r/an- 

nSZ 2 

This corresponds to two changes: Firstly, the original interaction vertex in ( 1431) allows 
all scattering processes such that 

2n 2 

Qri.si Qr 2 .s 2 ~l" Qr3,S3 Qr4,S4 H - k-i k/^ + IC4 G TL 

a 





n, 




r 2 ,s 2 


r 3 ,s 3 


r 4 ,s 4 


Restrictions 


1. 


r, 


s 


r, s 


r', s' 


r', s' 


(r,s)^(r',s'),s,s' = 0,±,2 


2. 


r, 


s 


r',s' 


r',s' 


r, s 


(r,s)^(7V),*,s' = 0,±,2 


3. 


r, 


s 


r, s 


r, s 


r, s 


s = 0,±,2 


4. 


r, 


s 


— r, s' 


— r, s 


r, s' 


( S ,s') = (±,T),(0,2),(2,0) 


5. 


r, 


s 


r, s' 


— r, s 


— r, s' 


( S , S ') = (±,T),(0,2),(2,0) 


6. 


r, 


s 


— r, s 


r, s 


— r, s 


s = 0,2 


7. 


r, 


s 


r',s' 


r, s 


-r',s' 


s = 0,2, s' = ± 


8. 


r, 


s 


r', s' 


— r, s 


r', s' 


s = ±, s ' = 0,2 


9. 


r, 


s 


r', s' 


r, s 


r', s' 


s,s' = 0,2 


10. 


r, 


s 


r', s' 


-r', s' 


— r, s 


(s,s') = (0,2),(2,0) 


11. 


r, 


s 


— r, s 


r', s' 


-r', s' 


( S , S ') = (0,2),(2,0) 



Table 1: List of all interactions terms potentially contributing to the effective model derived 
in Section |5T2| r, r' = ±. 
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but we now enforce, in addition, 

2tt 



Secondly, we modify the argument of the weight function u in (1431) by ignoring the depen- 
dence on ki — k2. Since the Q r , s are large and the (typically) small on the scale of 
our short-distance cutoff ir/a, we expect that the difference between the original- and our 
approximate interaction does not change the low energy properties of the model. 

We therefore only keep the interaction terms between the fermion flavors for which the 
condition in f|4*9|) is fulfilled. For Q = ir /2 we find 512 such terms, but for Q ^ 7r/2 there 
are only the 196 terms listed in table [TJ Below we discuss the corresponding interactions 
in more detail. 

Terms 1 in tabled] are Hartree-like, i.e. they describe scattering processes Q rs — > Q r ,s, 
Qr',s' — >■ Qr',s' which preserve fermion flavors. They are given by 

H l = E E(x) 6 «( )^ 1 -k2 + k3-k 4 ,0^(k 1 )^ s (k 2 )^, s ,(k 3 )i, s ,(k 4 ) 

r,s,r',s' kj 

where the prime in the sum indicates that the diagonal terms (r, s) = (r', s') are excluded 
(the latter are terms 3 in table [1] and are treated separately). It is useful to introduce the 
operators 

Mp)= E Cf) 2 ^U^)ks(^ 1+pM (50) 

ki,k2SA* s 

which allow to write the Hartree-like interactions in the following simpler form, 

H 'l = E E(f f V PrA-P)Pr'AP) (51) 



where we used (2n) 2 u(0) = a 2 V. The operators in fl50|) can be naturally interpreted as 
Fourier transformed fermion densities. 

Terms 2 are Fock-like, i.e. the corresponding scattering processes Q rjS — > Q r ', s '> Qr',s' — > 
Q, r s exchange fermion flavors. They are given by 

H 2 = E E(T) 6 «(Q^-Q^^ s 04 1 -k2 + k3-k4,o^ s (k 1 )i, s ,(k 2 )^ lS ,(k3)i lS (k 4 ). 

r,s,r',s' kj 

Renaming k 2 ,k 4 to k 4 , k 2 and using the anticommutator relations of the fermion field 
operators we can also write H' 2 in terms of fermion densities, 



#2 = ~ v r,s,r>, s> ^(z) 2 Pr A~~P) Pr>, s> (~P) + E V r,s,r' ,s' N r)S f r > jS > (52) 

r,s,r',s' p r,s,r',s' 

with N r>s = p r ,s(0) equal to the number operators in f l4"4"|) and the constants 

V r ,s,r',s' = (^) 2 u(Q r ,s ~ Qr' ,s') / ^ = ^r',s',r,s 



17 



and f rjS = XlkeA* i a /L) 2 . By simple computations we find 

u+,o,-,o = v +,2,-,2 = -V, u+,±,-,± = ^cos(2Q) 

V r ,0,r',± = Ur,0,r',2 = 0, ^r,±,r',=F = l^t 1 + COs(2Q)] (53) 
^r,±-,2 = VcOs(Q), V r ,±,+,2 = ~VcOs(Q), V r>s ^ s — V 

and 

fr,0 = ^ 2 , / r , ± = iK(l-K), / rj2 = I(l- K ) 2 (54) 

for all r, r' = ±. Note that / rjS equals the ratio of the area of A* s to the area of BZ, and 

Er,s fr,s = 1- 

We excluded the diagonal Hartree terms where (r, s) = (r', s') above since, by simple 
computations similar to the ones described above, they can be simplified to 

flS = £ E(x) 4 ^0)^Uk)^ s (k) = £ VN r , s f r , s (55) 

r,s k,k' r,s 

and thus do not contribute to the interactions; see Appendix [A] for details. 

Terms 4 and 5 in table Q] are mixed, i.e. they are Hartree-like in one and Fock-like in 
the other of the components of the momenta. By straightforward computations one finds 
that they exactly add up to zero; see Appendix [A] for details. 

Terms 6 are back-scattering terms where the condition in (149 j) holds true for a non-zero 
integer vector n. A simple computations shows that these terms are identical to zero; see 
Appendix lAl 

Terms 7-9 are BCS-like, i.e. Qri,si + Qr 3 ,s 3 £ 27rZ 2 /a and Q r2jS2 + Q r4jS4 € 2nZ 2 /a (note 
that Q rs + Q rs G 27rZ 2 /a for s = 0,2 and Q rjS + Q_ rjS = for s = ±). Straightforward 
computations show that all these terms are identically zero; see Appendix |A] There are 
also BCS-like terms (r 1 ,s 1 ) = (r, ±), (r 2 , s 2 ) = (-r, ±), (r 3 , s 3 ) = (-r, ±), (r 4 , s 4 ) = (r, ±) 
which give non-zero contributions, but they are already included in the Fock terms 2 above. 

Terms 10 and 11 are back-scattering terms. One can show that the former vanish, and 
the latter are equal to 

H 'n = ^EE(l) 4 (f ) 2 ^o(ki)^ r , (k 2 )^, 2 (k 3 )^ r , 2 (k 4 ); (56) 

kj r,r' 

see Appendix lAl 

For Q = 7r/2 the equation in (T49]) has many more solutions. This leads to additional 
interactions, including, for example, 

#bs = ^EEE(f) 4 (l) 2 i( k i)i r ,o(k2)*, s (k 3 )f r , s (k 4 ). 

kj r,r' s=± 

Nodal back-scattering terms like this, involving nodal fermion operators but which cannot 
be written in terms of nodal densities p r ,±{p), cannot be bosonized in a simple manner. In 
the following we therefore assume that Q ^ it/ 2. 
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To summarize, we found that, using the Approximations Al and A2 and assuming 
Q 7^ 7r/2, the free- and interaction parts of the 2D t-t'-V model Hamiltonian are equal to 

= E E(x) 2 [ e (Qns) + Mk) - ^J,.(k)^(k) (57) 

r,s k 

and H' = H[ + H' 2 + H' 3 + H' n , i.e. 

H' = E (V ~ V r,s,r>,s>) E(l) 2 Pn4-p)Pr', S '(p) + E V r,s,r>,s'N r , s f r , iS < + H' n (58) 
r,s : r f ,s' p r : s,r' : s' 

with the dispersion relations in (|47p and the constants in ([H]), (I46p . ( 15 3p and (15"4"|) . In the 
following we denote the full Hamiltonian H + if' as if + H'n (since H' u is not affected by 
the computations in the next section and will be dropped in Approximation A3'). 

Remark: The fact that many of the interaction terms listed in table [T] add up to zero is a 
consequence of the Pauli exclusion principle and holds true since we use spinless fermions. 

5.3 Normal ordering 

Up to now we assumed a reference state |0) in the fermion Fock space which is annihilated 
by all fermion operators ^(k). Before we can partly remove the UV cutoff below (Approx- 
imation A4) it is important to introduce another reference state Q in which all momentum 
states k inside some Fermi surface are filled. This state can be defined as follows, 

fi d = f n^(k)io) 

k 

where the product is over the set of all filled momenta k specified below. After Approx- 
imation A4 below one can only measure physical quantities with respect to this state Q, 
i.e., only expectation values of normal ordered operators are finite. We therefore have to 
determine the effect of normal ordering for various operators of interest to us. 

To be more specific, we assume that nodal fermion states (s = ±) are filled up to the 
Fermi surface arcs described in Section 13.21 (Hypothesis H2) and indicated by dashed lines 
in figure [TJ This is the only assumption we need to make about the Fermi surface since 
the partial continuum that we perform below only affects the nodal fermions. However, it 
is also convenient to assume that the antinodal fermions (s = and r = ±) have some 
filling v a /2 which is the same for r = + and r = — . We expect that there is an interesting 
parameter regime with v a = 1/2 [TU], but this remains to be shown [TT]. It is also convenient 
to assume that the in- and out fermions (s = 2 with r = — and r = +, respectively) are 
totally filled and totally empty, respectively. Thus the sets of filled states for the different 
fermion flavors r, s are as follows, 

S r ,± = {kG A; j± : rk± < 0} , S_, 2 = A*_ j2 , S + , 2 = 
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and SV,o such that the r.h.s. of (|33|) is the same for r = + and r = — (as will be seen, we 
need not be more specific about SV,o)- This allows us to fully characterize the new reference 
state as follows, 

$ )S (k)fi = VkeSV jS , Vv >s (k)fi = Vk£S r , s 

leading to fl32|) after Approximation A4 below. We compute the contributions v r>s = 
(a/L) 2 (Q, N r}S tt) of the different fermion flavors to the filling factor in this new reference 
state and find 

^0 = 1"/ Vt,± = i{1-k)(^-1 + k), u„ 2 = 1(1 - k) 2 , z/ +j2 = (59) 

with z/ a in ( l33i) . and for the total filling ^2 rs v r , s we obtain the result in (jHJ). 

We normal order various operators quadratic in the fermion fields as usual; see f l34|) . 
In particular, the normal ordered fermion number operators are :N riS : = N r ^ s — (L/a) 2 u ryS . 
Thus the filling is related to the expectation values of these normal ordered particle number 
operators as follows, 

" = 5 + (! -")£?- !) + *>- - \) + (f ) 2 E< :iV ^ : >■ ( 6 °) 

To compute the Hamiltonian in its normal ordered form we have to write it in terms 
of normal ordered fermion densities J r , s (p) = '■ Pr,s(p) : - Since normal ordering of these 
densities is irrelevant unless p = and p r , s {0) = N r s we obtain 

Jr,s(P) = PrAP) ~ S P,0 (f ) 2 Vr,s- (61) 

Inserting this into the Hamiltonian H obtained at the end of the previous section we get 
5=EE (¥) 2 [^(k)-^]:^(k)i iS (k): 

r,s keA* s 

' " (62) 

+ E ( V ~ V r,s,r',s>) E (f ) 2 ^(-P)^r', S '(p) + ^0 
r,s,r',s' p£ A* 

with the dispersion relations in (137]) . the parameters 

A*r,s = P — e(Qr,s) — 2V 2j ZV,s' — ^f r ,s,r',s'(/r',s' ~ 2lV,s')> (63) 

r',s' r',s' 

and the constants in ®, flU, §SD, flMD anc@ 

f ° = E E E (f) 2 KQr, s ) + ^(k)](fi,^ s (k)i )S (k)^}+^ (64) 

r=± s=0,±,2 keA* s 

9 The r, s summations in the following formulas are written explicitly since they are not changed in 
Approximation A3' below. 
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with 

S int = L 2 Vis 2 - \XV + ^ ^2 V r,s,r',s'Vr,s(fr',s' ~ U r> >') I (65) 
\ r,r'=± s,s'=0,±,2 / 

and i/ in dHJ). The constants /i rjS correspond to chemical potentials which are different for 
the different fermion flavors. 

The condition that the points Q r ,± are on the Fermi surface is equivalent to n T) ± = 0, 
and this fixes the chemical potential \x and \x r ^ s for r = ± and s = 0, 2. By straightforward 
computation we find /z± )0 = /xo i n an d 

M±,2 = =F[4t + V(l - «)(1 - k + 2(^2 - 1) cos(Q))][l ± cos(Q)] + At' sin 2 (Q). (66) 

The particle-hole transformation mentioned in Section H~T1 provides a useful check of our 
computations above: under this transformation, $ ra (k) -H- $- rg (k), and this is equivalent 
to the change of parameters in ( J23|) and 

Our approximate model above is invariant under this transformation since, in addition, 

Remark: A careful reader might wonder if one could also work with the normal ordered 
form of the interaction in (j2T|) . i.e. 

H' = J2 - y)^(x)^(yMyMx) 

(this would make the arguments in Section [5.21 somewhat simpler). The answer is "no" 
since Approximation A2 would spoil particle-hole transformation symmetry then. 

5.4 Partial continuum limit II 

We assume that the parameters are such that 

/i_ 2 < 0, /i +j2 > (67) 

on an energy scale of order t/5 or so. This gives some restrictions on the allowed parameter 
values that are important to check. If these conditions are fulfilled the in- and out fermions 
(s = 2 with r = — and r = +, respectively) have energies sufficiently far below and far 
above the Fermi energy, respectively, that the following approximation is expected to be 
appropriately 

A3': Drop all terms in the total Hamiltonian involving in- or out fermions. 

10 In the simplified treatment in [10] the following approximation was included in A4 below, and for this 
reason we label it as A3' rather than A3. 
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In particular we drop the interaction term H' n . Thus the total Hamiltonian is as in 
(|62|) . but from now on sums over r, s are restricted to r = ± and s = 0, ± unless stated 
otherwise. We thus obtain a model involving only nodal and antinodal fermions. 

We recall that J r ,s{p) is proportional to the sum of all terms : ipl s (ki)-^ riS (k 2 ) : with 
momenta kj e A* s such that p = k 2 — ki. Since the kj belong to finite sets the number of 
these terms gets smaller as \p±\ is increased, and it vanishes for p outside some finite set. 
In particular, J r , s =±(p) is non-zero only if |p s | < y/^nn/a and \p- s \ < v2(l —K)7r/a. In our 
Approximations A4 and A5 below we add terms to J r ,s=±(p)- The following approximation 
is to partly compensate for this. 

A3: Replace J r ,s=±(p) ^ n (EIP by Xs(p)Jr,s{p) with Xs(p) the cutoff functions in l[35\) . 

In the following approximation we partly remove the short-distance cutoff a for the 
nodal fermions by dropping the restriction on k s in A* s=± : 

A4: Replace the nodal Fourier spaces A* s=± in ( TjE ) by A* =± in l&U\). 

The final approximation adds terms to the nodal interactions and is necessary to obtain 
a model that can be bosonized in a simple manner. 

A5: Replace the nodal fermion densities thus obtained 

J r ,s=±(p) = Y,(t) 2 ^l(k-pNUk): (68) 
keA| 

by the ones in / fggj) . i.e., insert the integer sum. 

We thus obtain the Hamiltonian of the 2D Luttinger model as defined in Section 14.21 (To 
simplify our notation we introduced the constants in (jSJ) which we obtained by computing 
2a 2 (V — u r ,±,-r,±) = 9i etc.) Approximation A5 amounts to adding umklapp terms in the 
A;__ s -direction. Since these umklapp terms vanish in the full continuum limit a — > (at fixed 
parameters Vp etc.) we expect them to not (much) affect the low energy properties of the 
model. 

Remark: In Approximation A3 it is important to use some cutoff function Xs(p) to get 
a model that is mathematically well-defined. However, one has some freedom there: one 
could equally well replace the cutoff function in ( 135|) by 

\l if \Vs\< h -% and \Vs\<^^ 
I (J otherwise 

with additional free parameters bj > 0. According to Hypothesis HI in Section 13.21 the 
precise value of these parameters bj is not important. Our discussion above suggests bj < 2. 
Moreover, we find it natural to choose bi = b%, and the fact that J s =±,r(p) is defined for 
\P-s\ — (1 — k)tv / (\/2a) motivates our choice b 2 = 1. We note that the choice b\K = ^(l - K ) 
would be more convenient since it leads to simpler formulas. 
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5.5 Position space regularizations 



The fermion operators used in Section [2] are related to the ones in the present section by 
Fourier transformation 

^( x ) = ^E(x) 2 ^(k)e ik - x (69) 



2tt 
keA: 



The position vectors x therefore live on different spaces for different values of s: A* = ± is 
the Fourier space corresponding to spatial positions x where x s is continuous but X- s in a 
ID lattice Aid with lattice constant in (HI]) and L/a sites0 

Aid = {x = na : n G Z, |n| < L/(2a)} . 

The precise meaning of Hq = f d 2 x^2 r rvF :V4,±( X )(— i<9±)?/V,±(x) : is therefore 

nL/2 

Hq = dx± ^2 a^rvp :?/4 i± (x)(-i<9 ± )Vv,±( x ) : (70) 

J - L / 2 x T eA 1D r=± 

(Riemann sum in x T ) and the precise meaning of {ip r ±(x), V>±(y)} = 5 2 (x — y) is 

{Vv,±(x),V4,±(y)} = r&^^Oi -y±). (71) 

The reason for adding the umklapp terms in Approximation A5 is that (|69l) and (|39l for 
s = ± imply that 

J P , 8=± (x) = £ (J) 2 J r>s (p) e*' x =:< s (x)^ jS (x):, (72) 
pec s 

i.e. the nodal fermion densities are local. Note that there is another UV cutoff suppressed 
in the formal nodal interaction in (CQ), namely 

H' n = 2fd 2 xf d 2 y I 9l J2ks(* - y)J +>s (x) J_, s (y) + g 2 ^ 5+,_(x - y) J r . + (x) J r /,_(y) 

J \ a=± r,r'=± / 

with J d 2 x short for , 2 ctes Y1 X - eAi D ^ (^ or appropriate s) and 

Mx-y) = £ (i) 2 X,(p)x s '(p)e- ip - (x - y) 

peA* 

approximate Dirac delta functions. The formulas above allow to make the nodal Hamilto- 
nian in ([TJ precise as H n = + H~ +H' n . The precise meaning of the antinodal and mixed 
parts in (j2j) and can be given along the same lines but is not needed in the following. 



6 Bosonization and partial exact solution 

In this section we show that the 2D Luttinger model derived in the last section is amenable 
to analytical, non-perturbative computations. 

n Note that flUJ} ensures that L/(2a) £ (N + ±). 
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6.1 Bosonization 



The Hamiltonians in ( ITUj) can be interpreted as a model for ID fermions with x± a contin- 
uum spatial variable and x T a discrete flavor index. It is therefore possible bosonize the 
Hamiltonian using the same mathematical results that allow to bosonize the ID Luttinger 
model, as explained in more detail in Appendix |B] We collect these results, adapted to our 
needs, in the following. 

Proposition: (a) The nodal density operators are well-defined operators obeying the com- 
mutator relations 

[Jr,±{l>), Jr>,±(p')} = rSry—^-Sp^p'i^) 2 (73) 

and [J r ,±(p), Jr',^{p')] — 0. Moreover, J r ,±(p)^ = Jr,±(~P) an d 

J P>± (p)fi = Vp such that rp± > 0. (74) 

(b ) The following identity holds true 

£ (f frk ± :^ ±(k )^ ± ( k): = &7r £ (I) 2 : J r>± (-p) J r>± (p) : (75) 
kgA5_ pec± 

with C± in ( T77J[ ). 

(Proof outlined in Appendix IB1 with details in [12].) 

Part (a) of this proposition implies that the operators 

6 s= ±(±p) = Tt\ —— ^±,s(±p) Vp such that p s > 
are standard 2D boson operators, i.e. they obey the usual commutator relations 

etc. and b±(jp)Q = 0. Part (b) implies that the nodal Hamiltonian H n in ( )36|) is identical 
with a non-interacting boson Hamiltonian, and (|38|) shows that the antinodal densities are 
coupled linearly to these bosons. 

We note that the boson operators 6^2±(p) above are only defined if p s ^ 0, and H n 
includes also terms that cannot be expressed in terms of these boson operators due to the 
zero modes J r , s =±(p)\ Pa =o commuting with all boson operators. However, H n scales like L 2 
for large systems sizes L, whereas the sum of all these zero mode terms scales at most like 
L. This suggests that, for large systems, these zero mode terms can be ignored. To simplify 
our discussion below we therefore ignore zero mode terms, i.e. H n in the following differs 
from H n in fl36|) in that all terms with zero modes are implicitly excluded. A complete 
solution of the nodal model, including a detailed treatment of the zero modes, is given 
elsewhere [T2] . 

It is important to note that our scaling is such that the doping of the nodal fermions 
relative to Q is zero in the limit L/a — > oo; see fl97|) in Appendix [B] Moreover, the nodal 
filling v a is to be determined such that (a/L) 2 :N± t0 : = 0. Thus the doping constraint in 
(JSP]) can be simplified to the one in (jUJ). 
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6.2 Diagonalization of the nodal model 



A convenient way to write the bosonized nodal Hamiltonian is in terms of the operators 

iPs ^- 4n (Ps * 0) (76) 

n s (p) = y^[-i +lS (p) + i- )S (p)] 

which are standard 2D boson operators, i.e. they obey the usual canonical commutator 
relations 

[n s (p),^( P )] = -i^W(^) 2 

etc., with $s(p) = & s (~ p) an d nj,(p) = II S (— p). We refer to the quasi-particles corre- 
sponding to the fields in (lTB"j) as nodal bosons. 

Inserting this and f J75|) in f l36]) and fl38|) we obtain 

H » = f (E E (x) 2 : (t 1 - 7X.(p)]nl(p)ft.(p) + [1 + 7X.(p)]rf*J(p)*.(p)) = 

-± P ec s (7?) 

+ E (T) 2 7x + (P)X-(PW-S(p)i(p) + ^(P)^(P)] 

pec + nc_ 

and 

Hna = ^=Yl (i) 2 E^(-p)E^(p) 27ri ^^(p) ( 78 ) 

v pSC s r=± s=± 

with 7 = gi/ (vpair) = 2g 2 / '{vpcm). Inserting (jBJ) and (ITT]) we obtain (IT3j) . 

It is straightforward to diagonalize the nodal Hamiltonian in (ITT]) by a boson Bogoliubov 
transformation [12J. We obtain 

H n = £ n + E W.(P)S(P)6.(P) (79) 

a=± p£C s 

with 



^V^V ( |p|2 + VIP| 4 -^(2P+P-) 2 ) if \p±\ < nn min /(V2a) 
u s=± (p) = < WF |p s | if > 7TK/(V2a) (80) 

^-y/i — 7 2 |p s | otherwise 
the boson dispersion relation, where we use the shorthand notation 

def 7 2 def 
A = 1 ~ V9 ' = mill ( K > 1 - «)■ (81) 

(1 +7) 2 

The constant £ n is the nodal ground state energy, and we obtain 



2 

s=± pGC s 
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which is finite since oo s (p) ~ Vf\p s \ is non-zero only for a finite number of terms. The b^\p) 
are Bogoliubov transformed boson operators and can be computed explicitly [12]. From 
( 179]) is is easy to construct a complete set of exact eigenstates of H n . This solution shows 
that H n is well-defined if and only if 7 < 1, as mentioned earlier. 

Note that only a finite number of nodal boson modes are coupled to the antinodal 
fermions. Moreover, there are only a finite number of antinodal fermion degrees of freedom. 
Thus the 2D Luttinger Hamiltonian is well-defined if and only if H n is. 



6.3 Effective antinodal model 

In the previous section we found that the 2D Luttinger Hamiltonian can be mapped exactly 
to a Hamiltonian of non-interacting nodal bosons coupled to the antinodal fermions. This 
model is similar to the standard model of metallic electrons coupled to phonons. Similarly 
as for the latter model one can eliminate the nodal bosons and thus obtain an effective 
model for the antinodal fermions. 

We derive this effective model using a standard functional integral formalism |33j. To 
simplify our presentation in the main text we work in position space and suppress the UV 
and IR cutoffs as in Section [2j All formulas below can be made precise in Fourier space; 
see Appendix IC.ll 



We represent the partition function Z of the bosonized 2D Luttinger model in (jzj), (112]) 
and (fTB~]) as functional integral over the real- valued fields Q± = Q± (r, x) and Grassmann 
number valued fields ip± t0 = ip±fl(r, x) with r 6 [0, (3) the usual Matsubara time and (3 > 
the inverse temperature. Denoting the standard functional integral measures as D[- • ■] we 
obtain 

Z = Tr(e-P H ) = J D[<$> + , $_]D[V> + , 0> V'-.o] e^" 5 "- 5 ™ 
with the actions S(.) = L dr J d 2 x £(.) where 

^£ 1 



2 , ± 



v F (l - 7) 



(«9 T $ S ) 2 + M1 + 7 )(W 



\2 



+ v F y(d+$+)(d-$. 



with d T = -§p, 



£a = J2 $,0 ( d r + rc F d + d^ + c' F {d 2 + + d 2 _) - no) ^ + 2g 3 J+, J-,o (83) 



r=± 



with J ±)0 = Vi.oV'i.O' and 



=± 



By standard abuse of notation we denote the integration variables in the functional integral 
by the same symbol as the corresponding operators in the Hamilton formalism, and we 
suppress the common arguments r, x. 
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The boson integral is Gaussian and can be computed exactly, and the result has the 
following form, 

Z = Z n J D[^,o,V-,o]e- 5eff , Z n = J e~ s " (84) 

with — log Z n j (PL 2 ) the nodal contribution to the free energy and 

S c s = S a + J dr J d 2 x J dr' J d 2 y ^ Jo,r( r > x)w cfr (r - r, x - y) J y(r', y) (85) 



the effective action for the antinodal fermions, including the effective antinodal two-body 
potential v e g induced by the nodal bosons. We compute v e g and Z n in Appendix IC.U 

We obtain the following Fourier transformed effective potential^! 



v eS (uj,p) 



(g 4 ) 2 

Tiavp 



u 2 + u s (p) 2 



(86) 



with the functions 



u|(l-7) \p\ 2 + s 



(p 2 +-pD 2 +j\p\ 4 

' {l+7)y/\p\ 4 -A(2p+p-)< 



W s=± {p) = {2v 2 F {l- 1 )p 2 s 







if \p±\ < 7r« min /(\/2a) 

\p±\ < 7TK/(y / 2a) and 
\p±\ ^ 7iK min /(V2a) 



if 



otherwise 



(87) 

and w s =±(p) in flHUl) ; a; G 2ttZ/(3 are the usual boson Matsubara frequencies. Note that the 
second option in ( !87l) can only occur if k < 1/2 (since otherwise K m i n = k), and similarly 
in (EH below. 



As expected from the above-mentioned analogy with phonons, the effective interaction 
induced by the nodal bosons has a non-trivial time dependence. One can study the effec- 
tive antinodal fermion model defined by the action in (184"|) using a generalized mean field 
theory allowing for frequency dependent order parameters (this would be similar to the 
Migdal-Eliashberg theory of electron-phonon systems [3U [35]), but it is useful to have a 
simpler treatment. We therefore approximate the time dependent effective potential by an 
interaction local in time as follows, 



feff(T,X) 



6{t) I \ cS {r'^)dr' = 5(t) (i) 2 ^e ff (0,p)e i P- x 
Jo 



(88) 



peA* 



(arguments to justify this approximation can be found in the last paragraph of this section 
and in Appendix IC.2j) . In this time local approximation the effective antinodal action in 
flHU) can be obtained from the Hamiltonian 



H cS = H a + (l) 2 Yl M°.P) :X,o(-p)X',o(p) 
peA* 



39) 



r,r'=± 



12 i.e. v eS (r, X ) = l/(/3i 2 )^ p «( W ,p)e- i — j p- x 
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with H a in (137|) and the colons indicating normal ordering with respect to the state fl. 



We find that 



#eff(0,p)) 



-g eS if |p ± | < TTK min /(V2a) 

^ |b ± |<W(V2a) and 

otherwise 



with the constant g e s = (g^) 2 /[7iavp(l + 27)]. Inserting OH]) and (TTTT) we obtain c/ ef j- in 
(JTSJ). It is remarkable that v e s (0, p) is piecewise constant. 

Since J±,o(p) is non-zero only if \p±\ < 7r«/(\/2a), and -0 e fj(0, p) is constant and equal 
to — g e ff in that region for k < 1/2, one can simplify the effective Hamiltonian in ( 18 9 p for 
k < 1/2 as follows, 

H cS = H a - 2g eS (i) 2 ^+,o(-p)^-,o(p) (91) 
peA* 

(see Appendix IC.3I for details). We thus obtain the remarkable result that, for k < 1/2, 
the only effect of the nodal bosons in the time local approximation is a renormalization 
9s ^ 93 ~ 9cS of the antinodal coupling constant. 

We note that ttk min /(V2a) in (190]) can be regarded UV cutoff for the effective 
antinodal interaction, and according to Hypothesis HI in Section 13.21 we expect that it 
should be possible to remove this cutoff without significantly changing the low energy 
physics of the model. This suggests that we can simplify the effective antinodal model and 
use (]9~IT) even for n > 1/2. 

To justify the approximation in (188]) we note that it is similar to the one used by 
Bardeen and Pines to derive the phonon induced electron-electron interaction in a metal 
by eliminating the phonons in an electron-phonon Hamiltonian using a similarity transfor- 
mation [36J. Indeed, our bosonized 2D Luttinger Hamiltonian has the very same structure 
as an electron-phonon Hamiltonian, and one can therefore eliminate the nodal bosons in 
the Hamiltonian formalism using the method in [36]. This yields 

H?I = H a + (l) 2 (x) 4 E ^(k,p):^ (k)i )0 (k + p)^, ( k ')i, (k'- P ): 

p,k,k' r,r'=± 

with the following effective potential 

* n, \- M 2 \- W s (p) ( 1 



™f 2w s (p) \Us(p) - Ae r)0 (k, p) u 8 (p) + Ae r / )0 (k, p) 



depending on the energy differences Ae r (k, p) = £ r ,o(k) — £ ri o(k — p). However, different 
from electron-phonon systems, in the present case these energy differences are typically 
much smaller than w±(p) for |k|, |k'| <C ir/a (since £±.o(k) scales like ta 2 |k| 2 and w±(p) like 
to. I p I ) - Moreover, since W s (p) vanishes like (ta|p|) 2 for p — > 0, this effective interaction 
remains finite in this limit. This suggests that it is legitimate to approximate this interaction 
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by setting Ae: rj o(k, p) = 0. This yields our simplified time local approximation in (189]) . In 
the present case the fermion-boson interaction need not be small, however, and thus the 
argument in [36] is not conclusive. We therefore give a complementary argument to justify 
our approximation in Appendix IC.21 It would be interesting to adapt Migdal-Eliashberg 
theory [3U [35] to the present case and study the effective antinodal action without such 
approximation. 

7 Final remarks 

1. As mentioned in the introduction, the nodal Hamiltonian in ([T]) is essentially Mattis' 
model [21] . The relevance of the antinodal Hamiltonian in fl2]) as an effective model for 2D 
lattice fermion systems was emphasized by Schulz [25]. We argue in this paper that these 
two models capture complementary aspects of the physics of 2D lattice fermion systems, 
and our 2D Luttinger model includes them both. 

2. Together with de Woul we studied the effective antinodal Hamiltonian derived in Sec- 
tion 16.31 using mean field theory [11] . We found a significant parameter regime away from 
half-filling where the antinodal fermions are half filled [y a = 1/2) and fully gapped [ITj . 
We expect that the antinodal fermions do not affect low energy properties if they are fully 
gapped. We therefore propose that the exactly solvable Mattis model in (CQ) is an appro- 
priate low energy effective model in this regime. 

3. There exists also a parameter regime where the antinodal fermions are only partially 
gapped or not gapped at all [y a ^ 1/2) [11J. In this regime also the antinodal fermions 
contribute to the low energy properties of the system. We expect that this makes a qualita- 
tive difference in the physical behavior, and that the mean field treatment of the antinodal 
fermions is not as trustworthy for v a ^ 1/2 as for v a = 1/2. 

4. At very large doping values the condition in f )67|) eventually fails. In this regime we 
expect that one can ignore all but the in- or out fermions (s = 2 with r = — or r = +) for 
y < 1 or 1 - j/ < 1, respectively, and arguments similar to the ones in Section [5] lead to 
the following effective model for the 2D t-t'-V model, 

H= [(TCi,/2 + c' F )k 2 -/i ±i2 ]^ i2 (k)4 i2 (k) 

keA| 2 

(the interaction vanishes in our Approximation A2 due the the Pauli principle) . The renor- 
malized chemical potential ji2 t ± is determined by filling v. This effective model is trivially 
solvable. It suggests that the 2D t-t'-V model describes a Fermi liquid in the overdoped 
regime. 

5. The 2D Luttinger models has three parameters more than the 2D t-t'-V model, namely 
k, Q and v a . Our approach becomes more predictive if one can determine these parameters. 
As shown in [IT] , one can compute Q and v a self-consistently by using mean field theory 
for the antinodal fermions. This leaves only one free parameter, namely k. We expect that 
k can be fixed by comparison with experimental data or by the requirement that the mean 
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field phase diagram of the effective antinodal fermion Hamiltonian is consistent with the 
one of the 2D t-t'-V model. It would be interesting to fix k by a variational computation. 

6. Approximations A1-A5 in Section amount to adding or ignoring terms in the Hamil- 
tonian that, as we argued in Section I3.2[ are not relevant for the low energy properties of 
the system. It would be interesting to investigate in more detail if these approximations 
are relevant or not. This can be done, in principle, using our results: one can also bosonize 
the terms that were modified in our approximations (recall that the approximations that 
do not affect the nodal fermions can be easily avoided). It should be possible to investigate 
the relevance of these terms using renormalization group methods, similarly as in |37j . 

7. It is known that the truncated nodal model in (I3"oT) with g 2 = has "Luttinger liquid" 
behavior after |38J but not before Approximation A5 [2S]. We believe that our model with 
(72 > is less sensitive to Approximation A5 than the truncated model with g 2 = 0. The 
reason is that the boson dispersion relations in flHUj) are quite isotropic, whereas for g 2 = 
they are u>±(p) oc \p±\ independent of p T . We therefore expect that the boson propagator 
has better decaying properties in Fourier space for g 2 > than for g 2 = 0. 
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A Interaction terms 

This appendix contains details of the computation mentioned in Section 15.21 

To find all solutions of fH9|) we collected the vectors Q ri , Sl — Qr 2 ,s 2 in a 8 x 8 matrix 
with rows and columns labeled by (ri,Si) and (r 2 ,s 2 ), respectively. We used this matrix 
to find, for each case (ri, Si, r 2 , s 2 ), all cases (r 3 , S3, r 4 , s 4 ) such that f|49l) holds true. For 
Q = 71/ 2 we found 8 3 = 512 such cases. For Q ^ n/2 this number is reduced to the 196 
cases listed in table [TJ 

We now compute the interaction terms 3-11 in table [H We suppress sums over r, s (and 
r' and/or s' if applicable) when possible. Note that the possible values of r, s (etc.) are 
different for the different cases; see table [TJ 

The diagonal Hartree terms 3 are 

^3 = ED! )^(o)4 1 -k 2+k 3-k 4 ,o^(k 1 )^(k 2 )^t s (k 3 )^(k 4 ). 

r,s kj 

Rewriting this by renaming k 2 , k 4 to k 4 , k 2 , adding the result to H' 3 above, and dividing by 
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2 yields 

H * = EE(x) 6 «( )^- k2+ k3- k4 4^^ + 

r,s kj 

^ s (ki)^ r)S (k4)^ )S (k3)^ r)S (k 2 )J. 

We now use the canonical anticommutator relations of the fermion field operator twice and 
obtain the result in (I55|) . 

The sum of terms 4 and 5 can be written as 

#4+5 = E ii;Y U (Qr, S ~ Q_r, S ')^ki-k 2 +k 3 -k4,0 ^- r> s' ( k 2)^-r, S ( k 3) 4> r ,s' ( k *) 

+^, s (k 1 )i iS ,(k 2 )^l riS (k 3 )^_ riS ,(k 4 )) 

since u(Q r ,s — Q-r,s') = u(Q r>s — Q r , s ') for all cases (s,s') listed in table [Q Renaming in 
the second term k 2 , k 4 to k 4 , k 2 and using the canonical anticommutator relations of the 
fermion field operators twice we find H' 4+5 = 0. 

The back-scattering terms 6 are 

#6 = E (t) 6? HQ^ - Q-r,s)5 kl -k 2 +k3-k4,o ^ iS (ki)^_ riS (k 2 )^ r t s (k 3 )^_ riS (k4). 

Renaming k 1; k 3 to k 3 ,k! and using the canonical anticommutator relations we find that 
H' 6 = -H' 6 , i.e. H' 6 — 0. 

The BCS terms 7 and 9 are 

H i = E (t)^(Q^ - Qr^ s 0^k 1 -k 2+ k3-k 4 ,o^ s (k 1 )i /IS ,(k 2 )^t s (k 3 )^_ r , s/ (k 4 ) = o 

and 

^9 = E (X)^(Q^ - Qr',s05k 1 -k 2+ k 3 -k 4 ,O^ lS (ki)^ ) s ,(k 2 )^ s (k 3 )^ i s ,(k 4 ) = 

by the same argument as for Hq above. Since H' % is the Hilbert space adjoint of H' 7 this 
also proves H' 8 = 0. 

For the terms 10 we find u(Q r s — Qr' s') = for all cases listed in table [H and thus 
«io = °- 

Finally, the sum of all terms 11 is 

F n = EE ) 6 «(Q)^-k 2+ k3-k 4j o (^ i0 (k 1 )^ ri0 (k 2 )^, 2 (k 3 )^_ r , 2 (k 4 ) + 

kj r.r' 

^, 2 (ki)^_ r , 2 (k 2 )^, (k 3 )^_ r , (k 4 )). 
Using the canonical anticommutator relations and inserting (|28p we obtain (|56|) . 
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B Bosonization in ID 



We collect some well-known mathematical results that we need to bosonize the nodal 
fermion Hamiltonian. More details can be found in [2J [27J HQ], for example. 

Consider the Hamiltonian 

pL/2 M 

H = dxJ^Yl rV F'-Xt,n( X )(- [d )^r,n( X ) : ( 92 ) 

J-L/2 r=±n= _ M 

with Xr]n(x) ID fermion operators obeying canonical anti-commutator relations normalized 
such that 

{Xr,n{ x )i Xt',n'(y)} = Sn,n'$r,r>S(x ~ y), 

vp > a constant, d = d/dx, and n a flavor index; M > is some integer. We use 
anti-periodic boundary conditions, i.e. Fourier transform is 



f L/2 

Xr,n( k ) = (2tt)- 1/2 / dx X r,n( X ) ^ with ^|(Z + §) 
J-L/2 



and the colons indicate normal ordering with respect to a vacuum state Q defined by the 
following conditions, 

x r jm = xU-m = o \/rk>o. 

This Hamiltonian describes 2M + 1 flavors of ID Dirac fermions. We also introduce the 
fermion density operators 

jr,n{ X ) = ■Xt,n( X )Xr,n( X ) : • ( 93 ) 



We now state the well-known facts that allow to bosonize this Hamiltonian (see e.g. 
or [ID] for proofs): Firstly, the fermion densities obey the following non-trivial commutator 
relations, 

[jr,n( X )dr',n'(y)} = r5 r ,r'5 n . n > —05(x - y). (94) 

Z7T1 

Secondly, the Fourier transformed fermion densities 

3rM = J dxj r , n (x) e~ ipX , P G f (Z + |) 
Obey ] r ,n(py = Jr,n(-P) aIld 



>,„(p)ft = Vrp>0. (95) 

Finally, 

r 



L/2 pL/2 

dx :xt Jx)(-id)x r J x ) : = n / dx '■ JrA x ) 2 '■ ■ ( 96 ) 



-L/2 J-L/2 

We now observe that, by identifying x± with x, x T with nd, 2M + 1 = L/d, and 

Xr,n{x) = Va~il} r ,±(x), 
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the ID Luttinger Hamiltonian in (192]) is identical with the free nodal Hamiltonian in ( I70|) . 
We thus can identify 

and write ( l94l) as 

[J r ,±(x), J r ',±(y)] = roV >r /— r<9 ± o" 2 (x - y) 

Z7Tia 

with o~ 2 (x — y) = 8(x± — y±)S x - f) y- f /a. Moreover, summing the identity in fl96l) over n we 
get 

/ d 2 xr :ipl < ±(-x.)(—id±)ijj r ±('x)\= an J d 2 x :J r> ±(x) 2 : 

7j /2 

with J d 2 x = j_ L j 2 dx± J^x^eAiD ^' The results summarized in the proposition in Section 16. II 
are obtained from this and other identities above by Fourier transformation. 



We finally note that our scaling is such that 

w 

' — ( ' r) 

n=-M ' 



~ d f 1 M f 

j ■N r) ±:= I dx ± J ri± (x) = 2M ^ / dx j r<n (z 



are well-defined operators in the quantum field theory limit L/a — > oo, and therefore 



:t) 2 < 



:N rt ±:)^0 as L/a ^ oo. (97) 



C Effective interaction 

This section contains details on how the results presented in Section [6] were obtained. 
C.l Integrating out the nodal bosons 

To make precise the actions in the beginning of Section 16.31 we work in Fourier space. We 
suppress the arguments and write $± short for the Fourier transform Q±(u, p) of $±(r, x), 
with bj G 27rZ//3 the usual boson Matsubara frequencies. We also write K± short for 

K±(w,V>) = (™) _1/2 fi , 4ip±X±(p)L (c.;,p)/(27r) 
with J = J +> o + J-,o- Moreover, $±(u>, p) = <&±(—uj, — p), and similarly for K± and Jo- 
Using (177|) and (ITS]) we can give the following precise meaning to the action introduced 
at the beginning of Section 16.31 



with 



/3L 2 

w \s=±pec s pec 



w 2 

A s=± (w, p) = — — + + 7Xs(p))p 2 , B(u, p) = vp-fp+p- 

v F (l - 7XHPJ) 
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and 

C d = f (peA*: | P± |< 7rmin( ^ 1 - K) 

the region in Fourier space where the <f> + -<f>_ coupling term in (136 p is non-zero. 
It is convenient to introduce the matrix notation 



and similarly for i^. We also split the p-sums into sums over C and the complements 

cr={peC s : p^C}. 

This allows us to write 
2vr 2 



E ( E [it 1 ! + + + E E [^.*& + 2^I*J I 

w \pec s=± pe c ; J- / 



with 

,-i / A + B 



V 



B A. 



We can now compute the nodal boson functional integral by completing the square. This 
yields for the effective antinodal action defined in 



s eff = s a -^E(E + E E j-*J*. W«4EE (i) 2 ^ A 

« \ P eC s=± P ecJ- / w P eA* 

with 

f -(27rS)- 1 (( 7 4) 2 b + ,^)^,p)(p + ,P-) i if P e C 

v eS (u, p) = < -(27va)- 1 (g A ) 2 p 2 s /A s (uj, p) if p G Cj, s = + or - 

y otherwise 

and "t" indicating matrix transposition. This implies the result in fl86|) jff. We also obtain 



■ - mi ii . ; 



uj' z + oo s {p) 2 
lo s=± P ec 3 sy ^> 

with u; s= ±(p) in (I80j) . This expression for is formal but can be made well-defined by a 
multiplicative regularization as follows, 

7 = TTTTTT (2n7r//3 ) 2 + (r^) 2 tttt sinh 2 (A)i; F |p s |/2) 
* iUUi « 2 + ^P) 2 MpU sinh 2 (^(p)/2) W 

with /3 > arbitrary. This regularization amounts to multiplying Z n with an irrelevant 
(infinite) constant. The nodal free energy — log(Z n )/((3L 2 ) computed from this generalizes 
the nodal ground state energy £ n in (152"j) to finite temperature. 
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C.2 Locality of the effective interaction 



It is easy to see that v e g(r, x) decays like 1/B? for R = a/x 2 + (vft) 2 ^> a, different from 
and electron-electron interaction potential induced by phonons which decays like 1/R in 2D. 
This suggests that the former potential decays quickly in time and that (IBB!) is reasonable. 

To see in more detail how this potential v e g(T, x) behaves for L ^> |x| ^> a we compute 
it in the limits a — > and L — > oo. Using ( 1B61) . ( 1B7|) and ( IBD1) we obtain by inverse Fourier 
transformation 

v cff (r,x) = -ff eff 2 ] ,a /y I V 5 



^M 3 V JT^v f \t\ 



with the special function 



(l + 2 7 ) /- 27r rfx v- 



1 + 2 7 ) /^ dx v- [1 + ^(x)]Mx) 2 - 3x 2 cos 2 ( X - ?)] 



8tt(1 + T )(l - 7 2 ) 7o 2tt ^ [e s ( X ) 2 + x 2 cos 2 ( X - <p)] s 
where 



A = 1 — 7 2 /(l + 7) 2 , and = |x|(cos <p, sin <^). Plotting this function / 7 (</?, x) using 

MATLAB we find that it does not change much with 7, that it is positive for small values 
of x, and that it diverges like l/\x\ for x — > 0. Moreover, for x >~ 1 this function oszillates 
strongly with x and and is predominantly negative. Thus v e s(r, x) is sharply peaked and 
strongly attractive for |r| <C |x|/Vp, but for |r| > \x\/vp it is strongly oszillatory. The 
fact that J (1tv c s(t, x) = — g e ff<5 2 (x) shows that this effective interaction averages in time 
to zero everywhere except in x = 0. We thus expect that the local approximation in 
is reasonably accurate. 



C.3 Effective antinodal Hamiltonian 

This appendix contains details on the effective antinodal Hamiltonians in flB9|) and f l9Tj) . 
We first explain why we have to normal order in flB9l) with respect to Q, and then show 
how to obtain ([H]) from flBH . 

In our computation in Section 16.31 we use a functional integral formalism in which the 
antinodal fermions are represented by Grassmann numbers. It is important that these 
Grassmann numbers are defined using fermion coherent states with respect to the state Q 
(see e.g. [23]) since this allows us to do the computation in a way which is manifestly particle- 
hole symmetric. To convert the Hamiltonian to a functional integral we therefore normal 
order with respect to Q and then replace the fermion operators by Grassmann numbers 
dropping normal ordering. Similarly, to convert our result in ( 1B^|) to a Hamiltonian using 
the time local approximation in f )88p . we have to normal order the result with respect to 
the state Q. 
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For k < 1/2 the effective antinodal Hamiltonian in ( 189]) simplifies to 
H cS = H a - g cS ^ (l) 2 E : 4o(-p)^r',o(p) : • 

pGA* r,r'=± 

This differs from the result in ( l9Ti) by terms proportional to 

(*)± = E(l) 2 :4o(-p)4o( P ): 
p 

since : J±,o(p) : = ^±,o(p)- We now show that (*)± = by the Pauli exclusion principle. 
We compute (*)± above by inserting J±,o(p) = : P±,o(p) : an d (f50|) . We obtain 

(*)± = E (l) 2 (f ) 4 ^- k2+k 3-k 4 ,o :^,o(ki)Vi ±) o(k2)Vii ) o(k3)i 1± (k4): 

since :: A :: 5 :: = : AB :. We rewrite this by renaming k2, k4 to k4, k2. Under the normal 
ordering sign fermion creation and annihilation operators can be anti-commuted, and thus 
we find (*)± = — (*)±, i.e. (*)± = 0. This proves our result in fj91|) . 

The argument above also shows that fl9T|) is obtained for k > 1/2 provided that we 
approximate v e g by removing the UV cutoff, as discussed at the end of Section 16.31 
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